Synergism between remdesivir and ribavirin leads to SARS-CoV-2 extinction in cell culture.

BACKGROUND AND PURPOSE: There is a need for effective anti-COVID-19 treatments, mainly for individuals at risk of severe disease such as the elderly and the immunosuppressed. Drug repositioning has proved effective in identifying drugs that can find a new application for the control of coronavirus disease, in particular COVID-19. The purpose of the present study was to find synergistic antiviral combinations for COVID-19 based on lethal mutagenesis. EXPERIMENTAL APPROACH: The effect of combinations of remdesivir and ribavirin on the infectivity of SARS-CoV-2 in cell culture has been tested. Viral populations were monitored by ultra-deep sequencing, and the decrease of infectivity as a result of the treatment was measured. KEY RESULTS: Remdesivir and ribavirin exerted a synergistic inhibitory activity against SARS-CoV-2, quantified both by CompuSyn (Chou-Talalay method) and Synergy Finder (ZIP-score model). In serial passage experiments, virus extinction was readily achieved with remdesivir-ribavirin combinations at concentrations well below their cytotoxic 50 value, but not with the drugs used individually. Deep sequencing of treated viral populations showed that remdesivir, ribavirin, and their combinations evoked significant increases of the number of viral mutations and haplotypes, as well as modification of diversity indices that characterize viral quasi-species. CONCLUSION AND IMPLICATIONS: SARS-CoV-2 extinction can be achieved by synergistic combination treatments based on lethal mutagenesis. In addition, the results offer prospects of triple drug treatments for effective SARS-CoV-2 suppression.

SARS-CoV-2 mutant spectra as variant of concern nurseries: endless variation?

Introduction: SARS-CoV-2 isolates of a given clade may contain low frequency genomes that encode amino acids or deletions which are typical of a different clade. Methods: Here we use high resolution ultra-deep sequencing to analyze SARS-CoV-2 mutant spectra. Results: In 6 out of 11 SARS-CoV-2 isolates from COVID-19 patients, the mutant spectrum of the spike (S)-coding region included two or more amino acids or deletions, that correspond to discordant viral clades. A similar observation is reported for laboratory populations of SARS-CoV-2 USA-WA1/2020, following a cell culture infection in the presence of remdesivir, ribavirin or their combinations. Moreover, some of the clade-discordant genome residues are found in the same haplotype within an amplicon. Discussion: We evaluate possible interpretations of these findings, and reviewed precedents for rapid selection of genomes with multiple mutations in RNA viruses. These considerations suggest that intra-host evolution may be sufficient to generate minority sequences which are closely related to sequences typical of other clades. The results provide a model for the origin of variants of concern during epidemic spread─in particular Omicron lineages─that does not require prolonged infection, involvement of immunocompromised individuals, or participation of intermediate, non-human hosts.

Incipient functional SARS-CoV-2 diversification identified through neural network haplotype maps.

Since its introduction in the human population, SARS-CoV-2 has evolved into multiple clades, but the events in its intrahost diversification are not well understood. Here, we compare three-dimensional (3D) self-organized neural haplotype maps (SOMs) of SARS-CoV-2 from thirty individual nasopharyngeal diagnostic samples obtained within a 19-day interval in Madrid (Spain), at the time of transition between clades 19 and 20. SOMs have been trained with the haplotype repertoire present in the mutant spectra of the nsp12- and spike (S)-coding regions. Each SOM consisted of a dominant neuron (displaying the maximum frequency), surrounded by a low-frequency neuron cloud. The sequence of the master (dominant) neuron was either identical to that of the reference Wuhan-Hu-1 genome or differed from it at one nucleotide position. Six different deviant haplotype sequences were identified among the master neurons. Some of the substitutions in the neural clouds affected critical sites of the nsp12-nsp8-nsp7 polymerase complex and resulted in altered kinetics of RNA synthesis in an in vitro primer extension assay. Thus, the analysis has identified mutations that are relevant to modification of viral RNA synthesis, present in the mutant clouds of SARS-CoV-2 quasispecies. These mutations most likely occurred during intrahost diversification in several COVID-19 patients, during an initial stage of the pandemic, and within a brief time period.

Puzzles, challenges, and information reservoir of SARS-CoV-2 quasispecies.

Upon the emergence of SARS-CoV-2 in the human population, it was conjectured that for this coronavirus the dynamic intra-host heterogeneity typical of RNA viruses would be toned down. Nothing of this sort is observed. Here we review the main observations on the complexity and diverse composition of SARS-CoV-2 mutant spectra sampled from infected patients, within the framework of quasispecies dynamics. The analyses suggest that the information provided by myriads of genomic sequences within infected individuals may have a predictive value of the genomic sequences that acquire epidemiological relevance. Possibilities to reconcile the presence of broad mutant spectra in the large RNA coronavirus genome with its encoding a 3' to 5' exonuclease proofreading-repair activity are considered. Indeterminations in the behavior of individual viral genomes provide a benefit for the survival of the ensemble. We propose that this concept falls in the domain of "stochastic thinking," a notion that applies also to cellular processes, as a means for biological systems to face unexpected needs.

Fitness-Dependent, Mild Mutagenic Activity of Sofosbuvir for Hepatitis C Virus.

The concept of a mild mutagen was coined to describe a minor mutagenic activity exhibited by some nucleoside analogues that potentiated their efficacy as antiretroviral agents. In the present study, we report the mild mutagen activity of sofosbuvir (SOF) for hepatitis C virus (HCV). Serial passages of HCV in human hepatoma cells, in the presence of SOF at a concentration well below its cytotoxic concentration 50 (CC50) led to pre-extinction populations whose mutant spectra exhibited a significant increase of C→U transitions, relative to populations passaged in the absence of SOF. This was reflected in an increase in several diversity indices that were used to characterize viral quasispecies. The mild mutagenic activity of SOF was largely absent when it was tested with isogenic HCV populations that displayed high replicative fitness. Thus, SOF can act as a mild mutagen for HCV, depending on HCV fitness. Possible mechanisms by which the SOF mutagenic activity may contribute to its antiviral efficacy are discussed.

Hepatitis C virus fitness can influence the extent of infection-mediated epigenetic modifications in the host cells.

Introduction: Cellular epigenetic modifications occur in the course of viral infections. We previously documented that hepatitis C virus (HCV) infection of human hepatoma Huh-7.5 cells results in a core protein-mediated decrease of Aurora kinase B (AURKB) activity and phosphorylation of Serine 10 in histone H3 (H3Ser10ph) levels, with an affectation of inflammatory pathways. The possible role of HCV fitness in infection-derived cellular epigenetic modifications is not known. Methods: Here we approach this question using HCV populations that display a 2.3-fold increase in general fitness (infectious progeny production), and up to 45-fold increase of the exponential phase of intracellular viral growth rate, relative to the parental HCV population. Results: We show that infection resulted in a HCV fitness-dependent, average decrease of the levels of H3Ser10ph, AURKB, and histone H4 tri-methylated at Lysine 20 (H4K20m3) in the infected cell population. Remarkably, the decrease of H4K20m3, which is a hallmark of cellular transformation, was significant upon infection with high fitness HCV but not upon infection with basal fitness virus. Discussion: Here we propose two mechanisms ─which are not mutually exclusive─ to explain the effect of high viral fitness: an early advance in the number of infected cells, or larger number of replicating RNA molecules per cell. The implications of introducing HCV fitness as an influence in virus-host interactions, and for the course of liver disease, are warranted. Emphasis is made in the possibility that HCV-mediated hepatocellular carcinoma may be favoured by prolonged HCV infection of a human liver, a situation in which viral fitness is likely to increase.

Atypical Mutational Spectrum of SARS-CoV-2 Replicating in the Presence of Ribavirin.

We report that ribavirin exerts an inhibitory and mutagenic activity on SARS-CoV-2-infecting Vero cells, with a therapeutic index higher than 10. Deep sequencing analysis of the mutant spectrum of SARS-CoV-2 replicating in the absence or presence of ribavirin indicated an increase in the number of mutations, but not in deletions, and modification of diversity indices, expected from a mutagenic activity. Notably, the major mutation types enhanced by replication in the presence of ribavirin were A→G and U→C transitions, a pattern which is opposite to the dominance of G→A and C→U transitions previously described for most RNA viruses. Implications of the inhibitory activity of ribavirin, and the atypical mutational bias produced on SARS-CoV-2, for the search for synergistic anti-COVID-19 lethal mutagen combinations are discussed.

Viral Fitness, Population Complexity, Host Interactions, and Resistance to Antiviral Agents.

Fitness of viruses has become a standard parameter to quantify their adaptation to a biological environment. Fitness determinations for RNA viruses (and some highly variable DNA viruses) meet with several uncertainties. Of particular interest are those that arise from mutant spectrum complexity, absence of population equilibrium, and internal interactions among components of a mutant spectrum. Here, concepts, fitness measurements, limitations, and current views on experimental viral fitness landscapes are discussed. The effect of viral fitness on resistance to antiviral agents is covered in some detail since it constitutes a widespread problem in antiviral pharmacology, and a challenge for the design of effective antiviral treatments. Recent evidence with hepatitis C virus suggests the operation of mechanisms of antiviral resistance additional to the standard selection of drug-escape mutants. The possibility that high replicative fitness may be the driver of such alternative mechanisms is considered. New broad-spectrum antiviral designs that target viral fitness may curtail the impact of drug-escape mutants in treatment failures. We consider to what extent fitness-related concepts apply to coronaviruses and how they may affect strategies for COVID-19 prevention and treatment.

Efficacy decrease of antiviral agents when administered to ongoing hepatitis C virus infections in cell culture.

We report a quantification of the decrease of effectiveness of antiviral agents directed to hepatitis C virus, when the agents are added during an ongoing infection in cell culture vs. when they are added at the beginning of the infection. Major determinants of the decrease of inhibitory activity are the time post-infection of inhibitor administration and viral replicative fitness. The efficacy decrease has been documented with antiviral assays involving the combination of the direct-acting antiviral agents, daclatasvir and sofosbuvir, and with the combination of the lethal mutagens, favipiravir and ribavirin. The results suggest that strict antiviral effectiveness assays in preclinical trials may involve the use of high fitness viral populations and the delayed administration of the agents, relative to infection onset.

SARS-CoV-2 Mutant Spectra at Different Depth Levels Reveal an Overwhelming Abundance of Low Frequency Mutations.

Populations of RNA viruses are composed of complex and dynamic mixtures of variant genomes that are termed mutant spectra or mutant clouds. This applies also to SARS-CoV-2, and mutations that are detected at low frequency in an infected individual can be dominant (represented in the consensus sequence) in subsequent variants of interest or variants of concern. Here we briefly review the main conclusions of our work on mutant spectrum characterization of hepatitis C virus (HCV) and SARS-CoV-2 at the nucleotide and amino acid levels and address the following two new questions derived from previous results: (i) how is the SARS-CoV-2 mutant and deletion spectrum composition in diagnostic samples, when examined at progressively lower cut-off mutant frequency values in ultra-deep sequencing; (ii) how the frequency distribution of minority amino acid substitutions in SARS-CoV-2 compares with that of HCV sampled also from infected patients. The main conclusions are the following: (i) the number of different mutations found at low frequency in SARS-CoV-2 mutant spectra increases dramatically (50- to 100-fold) as the cut-off frequency for mutation detection is lowered from 0.5% to 0.1%, and (ii) that, contrary to HCV, SARS-CoV-2 mutant spectra exhibit a deficit of intermediate frequency amino acid substitutions. The possible origin and implications of mutant spectrum differences among RNA viruses are discussed.

SARS-CoV-2 Point Mutation and Deletion Spectra and Their Association with Different Disease Outcomes.

Mutant spectra of RNA viruses are important to understand viral pathogenesis and response to selective pressures. There is a need to characterize the complexity of mutant spectra in coronaviruses sampled from infected patients. In particular, the possible relationship between SARS-CoV-2 mutant spectrum complexity and disease associations has not been established. In the present study, we report an ultradeep sequencing (UDS) analysis of the mutant spectrum of amplicons from the nsp12 (polymerase)- and spike (S)-coding regions of 30 nasopharyngeal isolates (diagnostic samples) of SARS-CoV-2 of the first COVID-19 pandemic wave (Madrid, Spain, April 2020) classified according to the severity of ensuing COVID-19. Low-frequency mutations and deletions, counted relative to the consensus sequence of the corresponding isolate, were overwhelmingly abundant. We show that the average number of different point mutations, mutations per haplotype, and several diversity indices was significantly higher in SARS-CoV-2 isolated from patients who developed mild disease than in those associated with moderate or severe disease (exitus). No such bias was observed with RNA deletions. Location of amino acid substitutions in the three-dimensional structures of nsp12 (polymerase) and S suggest significant structural or functional effects. Thus, patients who develop mild symptoms may be a richer source of genetic variants of SARS-CoV-2 than patients with moderate or severe COVID-19. IMPORTANCE The study shows that mutant spectra of SARS-CoV-2 from diagnostic samples differ in point mutation abundance and complexity and that significantly larger values were observed in virus from patients who developed mild COVID-19 symptoms. Mutant spectrum complexity is not a uniform trait among isolates. The nature and location of low-frequency amino acid substitutions present in mutant spectra anticipate great potential for phenotypic diversification of SARS-CoV-2.

Vaccine breakthrough infections with SARS-CoV-2 Alpha mirror mutations in Delta Plus, Iota, and Omicron.

Replication of SARS-CoV-2 in the human population is defined by distributions of mutants that are present at different frequencies within the infected host and can be detected by ultra-deep sequencing techniques. In this study, we examined the SARS-CoV-2 mutant spectra of amplicons from the spike-coding (S-coding) region of 5 nasopharyngeal isolates derived from patients with vaccine breakthrough. Interestingly, all patients became infected with the Alpha variant, but amino acid substitutions that correspond to the Delta Plus, Iota, and Omicron variants were present in the mutant spectra of the resident virus. Deep sequencing analysis of SARS-CoV-2 from patients with vaccine breakthrough revealed a rich reservoir of mutant types and may also identify tolerated substitutions that can be represented in epidemiologically dominant variants.

In vitro Selection of High Affinity DNA and RNA Aptamers that Detect Hepatitis C Virus Core Protein of Genotypes 1 to 4 and Inhibit Virus Production in Cell Culture.

Hepatitis C virus (HCV) core is a highly conserved and multifunctional protein that forms the viral capsid, making it an attractive target for HCV detection and inhibition. Aptamers are in vitro selected, single-stranded nucleic acids (RNA or ssDNA) with growing applicability in viral diagnostics and therapy. We have carried out DNA and RNA in vitro selection against six different variants of HCV core protein: two versions of the full-length protein of genotype 1, and the hydrophilic domain of genotypes 1 to 4. The aptamer populations obtained were analyzed by means of Ultra-Deep Sequencing (UDS), the most abundant sequences were identified and a number of highly represented sequence motifs were unveiled. Affinity (measured as the dissociation constant, Kd) of the most abundant DNA and RNA aptamers were quantified using Enzyme-Linked OligoNucleotide Assay (ELONA)-based methods. Some aptamers with nanomolar or subnanomolar Kd values (as low as 0.4 nM) were the common outcome of DNA and RNA selections against different HCV core variants. They were tested in sandwich and competitive biosensor assays, reaching a limit of detection for HCV core of 2 pM. Additionally, the two most prevalent and high affinity aptamers were assayed in Huh-7.5 reporter cell lines infected with HCV, where they decreased both the viral progeny titer and the extracellular viral RNA level, while increasing the amount of intracellular viral RNA. Our results suggest that these aptamers inhibit HCV capsid assembly and virion formation, thus making them good candidate molecules for the design of novel therapeutic approaches for hepatitis C.

Guanosine inhibits hepatitis C virus replication and increases indel frequencies, associated with altered intracellular nucleotide pools.

In the course of experiments aimed at deciphering the inhibition mechanism of mycophenolic acid and ribavirin in hepatitis C virus (HCV) infection, we observed an inhibitory effect of the nucleoside guanosine (Gua). Here, we report that Gua, and not the other standard nucleosides, inhibits HCV replication in human hepatoma cells. Gua did not directly inhibit the in vitro polymerase activity of NS5B, but it modified the intracellular levels of nucleoside di- and tri-phosphates (NDPs and NTPs), leading to deficient HCV RNA replication and reduction of infectious progeny virus production. Changes in the concentrations of NTPs or NDPs modified NS5B RNA polymerase activity in vitro, in particular de novo RNA synthesis and template switching. Furthermore, the Gua-mediated changes were associated with a significant increase in the number of indels in viral RNA, which may account for the reduction of the specific infectivity of the viral progeny, suggesting the presence of defective genomes. Thus, a proper NTP:NDP balance appears to be critical to ensure HCV polymerase fidelity and minimal production of defective genomes.

Akt Phosphorylation of Hepatitis C Virus NS5B Regulates Polymerase Activity and Hepatitis C Virus Infection.

Hepatitis C virus (HCV) is a single-stranded RNA virus of positive polarity [ssRNA(+)] that replicates its genome through the activity of one of its proteins, called NS5B. This viral protein is responsible for copying the positive-polarity RNA genome into a negative-polarity RNA strand, which will be the template for new positive-polarity RNA genomes. The NS5B protein is phosphorylated by cellular kinases, including Akt. In this work, we have identified several amino acids of NS5B that are phosphorylated by Akt, with positions S27, T53, T267, and S282 giving the most robust results. Site-directed mutagenesis of these residues to mimic (Glu mutants) or prevent (Ala mutants) their phosphorylation resulted in a reduced NS5B in vitro RNA polymerase activity, except for the T267E mutant, the only non-conserved position of all those that are phosphorylated. In addition, in vitro transcribed RNAs derived from HCV complete infectious clones carrying mutations T53E/A and S282E/A were transfected in Huh-7.5 permissive cells, and supernatant viral titers were measured at 6 and 15 days post-transfection. No virus was rescued from the mutants except for T53A at 15 days post-transfection whose viral titer was statistically lower as compared to the wild type. Therefore, phosphorylation of NS5B by cellular kinases is a mechanism of viral polymerase inactivation. Whether this inactivation is a consequence of interaction with cellular kinases or a way to generate inactive NS5B that may have other functions are questions that need further experimental work.

A Two-Level, Intramutant Spectrum Haplotype Profile of Hepatitis C Virus Revealed by Self-Organized Maps.

RNA viruses replicate as complex mutant spectra termed viral quasispecies. The frequency of each individual genome in a mutant spectrum depends on its rate of generation and its relative fitness in the replicating population ensemble. The advent of deep sequencing methodologies allows for the first-time quantification of haplotype abundances within mutant spectra. There is no information on the haplotype profile of the resident genomes and how the landscape evolves when a virus replicates in a controlled cell culture environment. Here, we report the construction of intramutant spectrum haplotype landscapes of three amplicons of the NS5A-NS5B coding region of hepatitis C virus (HCV). Two-dimensional (2D) neural networks were constructed for 44 related HCV populations derived from a common clonal ancestor that was passaged up to 210 times in human hepatoma Huh-7.5 cells in the absence of external selective pressures. The haplotype profiles consisted of an extended dense basal platform, from which a lower number of protruding higher peaks emerged. As HCV increased its adaptation to the cells, the number of haplotype peaks within each mutant spectrum expanded, and their distribution shifted in the 2D network. The results show that extensive HCV replication in a monotonous cell culture environment does not limit HCV exploration of sequence space through haplotype peak movements. The landscapes reflect dynamic variation in the intramutant spectrum haplotype profile and may serve as a reference to interpret the modifications produced by external selective pressures or to compare with the landscapes of mutant spectra in complex in vivo environments. IMPORTANCE The study provides for the first time the haplotype profile and its variation in the course of virus adaptation to a cell culture environment in the absence of external selective constraints. The deep sequencing-based self-organized maps document a two-layer haplotype distribution with an ample basal platform and a lower number of protruding peaks. The results suggest an inferred intramutant spectrum fitness landscape structure that offers potential benefits for virus resilience to mutational inputs.

Historical Perspective on the Discovery of the Quasispecies Concept.

Viral quasispecies are dynamic distributions of nonidentical but closely related mutant and recombinant viral genomes subjected to a continuous process of genetic variation, competition, and selection that may act as a unit of selection. The quasispecies concept owes its theoretical origins to a model for the origin of life as a collection of mutant RNA replicators. Independently, experimental evidence for the quasispecies concept was obtained from sampling of bacteriophage clones, which revealed that the viral populations consisted of many mutant genomes whose frequency varied with time of replication. Similar findings were made in animal and plant RNA viruses. Quasispecies became a theoretical framework to understand viral population dynamics and adaptability. The evidence came at a time when mutations were considered rare events in genetics, a perception that was to change dramatically in subsequent decades. Indeed, viral quasispecies was the conceptual forefront of a remarkable degree of biological diversity, now evident for cell populations and organisms, not only for viruses. Quasispecies dynamics unveiled complexities in the behavior of viral populations,with consequences for disease mechanisms and control strategies. This review addresses the origin of the quasispecies concept, its major implications on both viral evolution and antiviral strategies, and current and future prospects.

Mutation Rates, Mutation Frequencies, and Proofreading-Repair Activities in RNA Virus Genetics.

The error rate displayed during template copying to produce viral RNA progeny is a biologically relevant parameter of the replication complexes of viruses. It has consequences for virus-host interactions, and it represents the first step in the diversification of viruses in nature. Measurements during infections and with purified viral polymerases indicate that mutation rates for RNA viruses are in the range of 10-3 to 10-6 copying errors per nucleotide incorporated into the nascent RNA product. Although viruses are thought to exploit high error rates for adaptation to changing environments, some of them possess misincorporation correcting activities. One of them is a proofreading-repair 3' to 5' exonuclease present in coronaviruses that may decrease the error rate during replication. Here we review experimental evidence and models of information maintenance that explain why elevated mutation rates have been preserved during the evolution of RNA (and some DNA) viruses. The models also offer an interpretation of why error correction mechanisms have evolved to maintain the stability of genetic information carried out by large viral RNA genomes such as the coronaviruses.

Population Disequilibrium as Promoter of Adaptive Explorations in Hepatitis C Virus.

Replication of RNA viruses is characterized by exploration of sequence space which facilitates their adaptation to changing environments. It is generally accepted that such exploration takes place mainly in response to positive selection, and that further diversification is boosted by modifications of virus population size, particularly bottleneck events. Our recent results with hepatitis C virus (HCV) have shown that the expansion in sequence space of a viral clone continues despite prolonged replication in a stable cell culture environment. Diagnosis of the expansion was based on the quantification of diversity indices, the occurrence of intra-population mutational waves (variations in mutant frequencies), and greater individual residue variations in mutant spectra than those anticipated from sequence alignments in data banks. In the present report, we review our previous results, and show additionally that mutational waves in amplicons from the NS5A-NS5B-coding region are equally prominent during HCV passage in the absence or presence of the mutagenic nucleotide analogues favipiravir or ribavirin. In addition, by extending our previous analysis to amplicons of the NS3- and NS5A-coding region, we provide further evidence of the incongruence between amino acid conservation scores in mutant spectra from infected patients and in the Los Alamos National Laboratory HCV data banks. We hypothesize that these observations have as a common origin a permanent state of HCV population disequilibrium even upon extensive viral replication in the absence of external selective constraints or changes in population size. Such a persistent disequilibrium-revealed by the changing composition of the mutant spectrum-may facilitate finding alternative mutational pathways for HCV antiviral resistance. The possible significance of our model for other genetically variable viruses is discussed.